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In nuclear collisions at RHIC energies, an excess of Q hyperons over 2 is observed, indicating 
that Q carries a net baryon number despite s and § quarks being produced in pairs. The baryon 
number in 2 could have been transported from the incident nuclei and/or acquired and balanced 
in baryon pair productions associated with other types of anti-hyperons, such as =. To investigate 
these two scenarios, we propose to measure correlations between Q and K, as well as between Q and 
anti-hyperons. We will use two versions, the default and string-melting, of a multiphase transport 
(AMPT) model to illustrate the correlation method. We will present the Q-hadron correlations from 
simulated Au+Au collisions at VJSNN = 7.7 and 14.6 GeV, and discuss the dependence on collision 
energy and on the hadronization scheme in these two AMPT versions. These correlations from the 
AMPT model provide a baseline for experimental exploration of the dynamics of baryon number 
transport and the effects of baryon number and strangeness conservation in nuclear collisions. 


I. INTRODUCTION 


Strangeness enhancement was proposed as a signature 
of the quark-gluon plasma (QGP) created in relativistic 
heavy-ion collisions [1] and has been a subject of intensive 
theoretical and experimental investigations [2]. Since the 
incident protons and neutrons are composed of u and d 
quarks, the strange quarks observed in the aftermath of 
the collision can only originate from pair of ss produc- 
tion. Lattice Quantum ChromoDynamics (QCD) pre- 
dicted that the temperature at which the quark-hadron 
phase transition occurs is approximately 150 MeV. Thus, 
in the QGP phase where the temperature is higher than 
the s-quark mass, strangeness may be abundantly pro- 
duced via flavor creation (qq > s5, gg — s5) and 
gluon splitting (g > s5), leading to enhanced production 
of strangeness in the final state [1, 3]. The investiga- 
tion of multi-strange hyperon production is particularly 
valuable for studying the equilibration of strangeness in 
the QGP. The yields of Q hyperons in heavy-ion colli- 
sions, for example, have been measured to be significantly 
higher than those in p+p collisions scaled by the num- 
ber of participants at CERN Super Proton Synchrotron 
(SPS) [4], Relativistic Heavy Ion Collider (RHIC) [5, 6], 
and the Large Hadron Collider (LHC) [7, 8]. Further- 
more, many strange hadrons are believed to have small 
hadronic rescattering cross-sections, so that they retain 
information from the hadronization stage and can be 
used to probe the phase boundary of the quark-hadron 
transition [9-11]. 

The RHIC Beam Energy Scan (BES) program aims 
to search for a possible critical point in the QCD phase 
diagram where strangeness production also plays a ma- 
jor role [12]. At low BES energies, the measured ratios 
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of anti-baryons to baryons are significantly lower than 
unity at midrapidities for three hyperons (A, =, and Q). 
Therefore, these hyperons must carry a net baryon num- 
ber. The baryon number transport dynamics have been 
a subject of interest of the heavy-ion collision physics 
since its inception [13-15]. The baryon number in a pro- 
ton or neutron may be attributed to valence u and d 
quarks, each carrying 1/3 of the baryon number. How- 
ever, in high energy collisions the valence quarks tend to 
inherit a significant fraction of the incident nucleon mo- 
mentum, making them ineffective in transporting baryon 
number from beam rapidities to midrapidities. Previous 
theoretical calculations have proposed the existence of 
topological objects known as gluon junctions, which can 
effectively convey baryon number over large rapidity gaps 
in nuclear collisions [16]. Once a gluon junction reaches 
midrapidities during a nuclear collision, it must emerge 
as a baryon in the final state, with its flavor determined 
by the quark flavors present in the surrounding QGP 
medium. Such exotic dynamics have a particularly dis- 
cernible impact on net 2 hyperons, as 2 hyperons consist 
of three s quarks that must be pair-produced. 

The quantum numbers of strangeness and baryon num- 
ber are strictly conserved in nuclear collisions, leading to 
correlations among particles in the final state. To gain 
insights into the production dynamics of Q hyperons, we 
propose to measure the correlations between Q and other 
particles, namely K+, A, and St, along with their re- 
spective anti-particle pairs. By examining the shape and 
strength of these correlation functions, we illustrate the 
extent that we can quantitatively characterize the role of 
conservation laws in the Q production dynamics in nu- 
clear collisions. These measurements will also provide a 
baseline reference to the search for exotic dynamics such 
as gluon junction transport. 

In this paper, we use a multiphase transport (AMPT) 


model [17] to simulate Au+Au collisions at \/syn = 7.7 
GeV and 14.6 GeV and present the Q-hadron correla- 
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FIG. 1. Schematic illustration of two possible scenarios for the Q7 production based on the quark coalescence picture with 


strangeness and baryon number conservation. In Scenario 1 (left), Q7 


carries the baryon number initially residing in the u and 


d quarks from colliding nuclei. In Scenario 2 (right), Q7 does not carry a net baryon number. Associated production of Q and 
kaon in default version could have similar qualitative features as scenario 1 in coalescence picture. 


tions from these simulations. In section 2, we describe 
briefly the AMPT simulations and the analysis methods. 
The correlation results and discussions are presented in 
section 3, followed by a summary in section 4. 


II. AMPT SIMULATIONS AND ANALYSIS 
METHODS 


A. AMPT Simulations 


To investigate the impact of different hadronization 
schemes on baryon number transport dynamics, we ex- 
ploit both the default and string melting (SM) versions 
of the AMPT model to simulate Au+Au collisions. The 
initial phase space is provided by the HIJING model [18- 
21], which is a Monte Carlo event generator for parton 
and particle production in high-energy hadronic and nu- 
clear collisions. In the default version of AMPT, minijets 
and their parent nucleons form excited strings after par- 
tonic interactions of minijets, and these strings fragment 
into hadrons via Lund String Fragmentation [22]. In the 
SM version, the strings are converted through the ” string 
melting” mechanism into partons, which subsequently 
interact with each other during the evolution, and the 
coalescence mechanism is used to combine partons into 
hadrons. The parton interactions are described by the 
Boltzmann equation, and solved by the ZPC model [23], 
which only includes two-body elastic scatterings. In the 
quark coalescence process, the two or three nearest par- 
tons (quarks and anti-quarks) in the phase space recom- 
bine to form a meson or a baryon. The AMPT version 
vl.25t4cu/v2.25t7cu, which strictly conserves net electric 
charge, strangeness, and baryon number, is employed in 
our study. We have generated more than 50 million min- 
imum bias events of Au+Au collisions at both 7.7 GeV 
and 14.6 GeV. 


B. Q Production and Conservation of Strangeness 
and Baryon Number 


The Q` production is governed by both strangeness 
conservation (SC) and baryon number conservation. Ad- 
ditionally, the presence of net Q hyperons at midrapidi- 
ties indicates the involvement of baryon number trans- 
port (BNT) dynamics. Figure 1 illustrates two basic sce- 
narios based on the quark coalescence picture for the Q7 
production. In Scenario 1 (left side of Fig. 1), three 5 
quarks are pair-produced with the three s quarks in Q7, 
and may combine with u or d quarks from the incident 
nuclei to form three kaons. In this case, the Q7~ hyperon 
does not have an accompanying anti-baryon and carries 
a baryon number initially residing in the u and d quarks 
from the colliding nuclei. In exotic dynamics involving 
gluon junction, the production of the three s-5 pairs may 
be combined with a gluon junction, which gives rise to 
Q-. Consequently, Scenario 1 would describe a BNT 
via gluon junction to the Q hyperon, while valence u/d 
quarks form the kaons [15]. In default version, associated 
production of Q and kaons could also yield qualitatively 
similar features as the coalescence picture. Therefore, 
Scenario 1 encompasses contributions from both ordi- 
nary physical processes and the gluon junction dynamics. 
The AMPT model in our simulation does not invoke the 
gluon junction dynamics, and our results may serve as a 
baseline to the experimental search for the exotic dynam- 
ics. In Scenario 2 (right side of Fig. 1), the three pair- 
produced 3 quarks combine with u(ū) or d(d) quarks to 
form an anti-hyperon and a kaon. For instance, =+(=°) 
and K°(K*) could emerge alongside the Q~ hyperon. In 
this scenario, the baryon number in Q` is balanced by 
the other types of anti-hyperon, and no direct BNT from 
the incident nuclei is present. 

To characterize the numbers of kaons and anti-baryons 
associated with the Q production, we introduce 


ANa = (A)w.9- = (A)w.o.9-, (1) 


where (A)w.o- and (A)w.o.n- denote the average num- 
bers of particle A in events with one Q7 and without any 


Q, respectively. Here, we assume that all other aspects 
of the two event classes are the same. Table I lists the 
ANx and ANg values expected by the two scenarios, 
where K represents both K+ and K®, and B refers to 
anti-baryons (A, ©, =, and so on) associated with the 
Q- production. In our study, we choose the beam en- 
ergies of 7.7 and 14.6 GeV in order to limit the fraction 
of Au+Au events that produce multiple Q (Q) hyperons. 
Compared with Scenario 2, Scenario 1 exhibits a stronger 


Q--K correlation and a weaker 2~-B correlation. 


TABLE I. ANx and ANg in the two scenarios for the pro- 
duction of an Q`. 


ANx AN5 
Scenario 1 (SC+BNT) 3 0 
Scenario 2 (SC) 1,2,3 1 


C. Correlation Using the Event Mixing Technique 


To divide out the combinatorial background, we first 
analyze the correlations between Q7 and strange hadrons 
with the traditional normalization using mixed events. 
The correlation function is 


B(k*)’ 2) 


where A(k*) is the same-event distribution, B(k*) is the 
mixed-event distribution, and k* = lki — kə| is the re- 
duced momentum in the pair rest frame. M is the nor- 
malization factor determined by matching the same-event 
and mixed-event correlations in the uncorrelated phase 
space, for example, by requiring C(k* > 1 GeV/c) = 1. 
The technique of event mixing normalization allows for 
the investigation of the correlation length between two 
types of particles in momentum space by analyzing their 
distributions. The underlying assumption is that the dis- 
tribution in mixed events should be equivalent to that in 
same events in the absence of any correlation, thereby 
representing the combinatorial background. By dividing 
out this mixed-event background, one can obtain the cor- 
relation between the two particles of interest. It is crucial 
to choose a proper k* range for normalization in order to 
achieve meaningful results. Typically, the k* range with 
the high counting density is selected for normalization, as 
it provides a reliable estimate of the background. This 
method has also been applied to examine the correla- 
tion function for Q~-K+, Q--A°, and Q--S+ in terms 
of relative rapidity and relative transverse momentum. 
However, it is important to note that due to the normal- 
ization procedure, this correlation function cannot fully 
capture the possibility of multiple kaons being correlated 
with the 2 production, nor can it fully account for the 
sensitivity of the correlation to the production dynamics. 


D. Correlations Using Combinatorial Background 
Subtraction 


In order to make the correlation measurement sensi- 
tive to the number of associated hadrons, we introduce 
the combinatorial background subtracted (CBS) corre- 
lations, by taking the difference between the Q-hadron 
correlation and the Q-hadron correlation, each normal- 
ized by the corresponding number of Q or Q hyperons. 
For example, the Q~-K* correlation is defined as 


= dNo-x+/dk*  dNga+g+/dk* 
o No Na ` 


Ca ko (k*) (3) 
This background subtraction approach is intended to ex- 
tract the main component of the correlation due to SC 
in the Q--Kt pairs. The opposite-sign pair (Q7-K+ 
or Q*+-K~—) distribution in the first term contains the 
signal, while the same-sign pair distribution in the sec- 
ond term (Q+-K+ or QT-K7) models the uncorrelated 
background. This subtraction scheme is sensitive to the 
difference in the number of kaons between events with 
QT and events with Qt, as well as to the phase space 
distribution of the extra kaons. There may be variations 
of kinematic phase spaces for these pairs in nuclear col- 
lisions, and such effects can be mitigated with this nor- 
malization scheme and with fine collision centrality bins. 
Similar approaches have been applied to study the cor- 
relations between Q and other hadrons. 


II. CORRELATION RESULTS AND 
DISCUSSIONS 


A. Strangeness Conservation and Strange Hadron 
Yields 


We first list in Table II the AMPT simulations of the 
difference in the number of s5 pairs between events with 
one QT and events without any Q7 or QF in Au+Au 
collisions at ysyyN = 7.7 GeV and 14.6 GeV. For events 
with one Q, at least three s5 pairs are produced because 
of SC. The AMPT results are slightly greater than three, 
and the excess indicates that the underlying strangeness 
production may affect the probability for the Q forma- 
tion. Compared with the SM version of AMPT, the de- 
fault version yields a slightly larger number of s5 pairs, as 
well as a smaller 2 formation probability. This correla- 
tion is presumably due to the difference in the formation 
dynamics and/or the s5 phase space. The AMPT simula- 
tions suggest that at the lower collision energies, slightly 
more strangeness in the underlying event is needed to 
form Q hyperons than at higher energies. 

Table III shows the AMPT calculations of ANg and 
AN g as defined in Eq. (1). When compared with the 
expectations in Table I corresponding to the two Q pro- 
duction scenarios, these numbers indicate that the Q pro- 
duction is likely to receive contributions from both sce- 
narios. The SM version seems to favor the dominance of 
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FIG. 2. AMPT calculations of the difference in strange hadron yields between Q7 events 


collisions at 7.7 GeV. 
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FIG. 3. AMPT calculations of the difference in strange hadron yields between Q7 events 


collisions at 14.6 GeV. 


TABLE II. Difference in the number of s5 pairs between 
events with one QT and events without any Q7 or QF. 


7.7 GeV 14.6 GeV 


3.22 3.07 
3.25 3.16 


AMPT SM 
AMPT default 


TABLE III. Average numbers of ANx and AN g as defined 
in Eq. (1). 


7.7 GeV 14.6 GeV 

ANx AN; ANx AN; 

AMPTSM 2.46 0.017 2.44 0.119 
AMPT default 1.74 0.078 1.76 0.28 


Scenario 1, with the ANg values close to three and the 
lower AN g values. 

Strangeness and baryon number can be represented by 
various particle types in the final states of nuclear colli- 
sions. For example, s quarks can exist in K+ and K®. 
The actual distributions of strangeness and baryon num- 
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bers in the final-state particles could be sensitive to nu- 
clear dynamics and may also depend on beam energy. 
Figures 2 and 3 show the AMPT results of the difference 
in strange hadron yields between Q7 events and non-Q— 
events of Au+Au collisions at 7.7 GeV and 14.6 GeV, 
respectively. At 7.7 GeV and 14.6 GeV, the correlations 
between Q7 and kaons are stronger in the SM version, 
but the correlations between 27 and anti-hyperons are 
stronger in the default version. For Q~ events, the num- 
ber of kaons is much larger than the number of anti- 
hyperons, suggesting that the correlations between Q7 
and kaons are much stronger than those between Q7 and 
anti-hyperons. These numbers seem to support that sce- 
nario 1 contributes more significantly to Q yields in both 
versions. 


B. Correlations between Q= and Strange Hadrons 


We first show the correlations between Q and strange 
hadrons as a function of the reduced momentum k* in 
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FIG. 4. Correlation functions between Q7 and strange hadrons using the event mixing normalization in the 0-5% centrality 


range of Au + Au collisions at 7.7 GeV (|n| < 1). 
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FIG. 5. Correlation functions between Q7 and strange hadrons using the event mixing normalization in the 0-5% centrality 


range of Au + Au collisions at 14.6 GeV (|n| < 1). 


the pair rest frame using the event mixing normaliza- 
tion. Figure 4 presents the correlation functions for (a) 
Q--K+, (b) Q--A, and (c) Q--E* pairs for 0-5% cen- 
trality Au + Au collisions at 7.7 GeV. The normalization 
between the same event and the mixed event distribu- 
tions is determined by the k* region of 0.6-1.5 GeV/c. 
There is a distinct correlation between Q7 and K*, in- 
dicating that SC must play a major role in the Q7 and 
K* yields, with the typical correlation length on the or- 
der of less than k* of 0.5 GeV/c. The QT may also be 
correlated with the anti-hyperons of A and Et, but the 
current simulated results do not have sufficient statistics 
to be definitive. 


Figure 5 shows (a) QT-K+, (b) Q--A, and (c) Q7-E* 
correlations for 0-5% centrality Au + Au collisions at 
14.6 GeV. The Q7-K* correlation at 14.6 GeV seems to 
be less prominent than that at 7.7 GeV. This could be 
explained by our normalization scheme that involves the 
dilution of the Q~-K* correlation by uncorrelated kaons, 
which are expected to be more significant at 14.6 GeV. 


According to Table I, the correlation between Q and 
anti-hyperons may also be susceptible to the two Q pro- 
duction scenarios. The middle and right panels of Figs. 4 
and 5 show the Q-A? and Q--S* correlation functions 
using the event-mixing technique at 7.7 GeV and 14.6 
GeV, respectively. At both energies, the Q~-A° results 


lack enough statistics for a definitive conclusion, whereas 
some level of the Q~-E* correlation may exist. However, 
there is no significant difference between Q~-E+ and Qt- 
~, similar to the 0-K results. 
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For both beam energies, there is no significant differ- 
ence in the observed correlations between the two AMPT 
hadronization schemes. It seems that the correlations us- 
ing the event-mixing technique are only sensitive to the 
kinematic region where SC links the ss pairs. The mag- 
nitude differences as shown in Table III are, however, not 
quantitatively reflected in the measured correlations due 
to the normalization scheme. We will explore the CBS 
correlations within the same event framework and back- 
ground subtraction scheme. 


Our goal is to use the Q+-K~ correlation (Scenario 2 
only) as a baseline for the Q~-K* correlation. For each 
of these CBS correlations, we again select events with 
one QT or Q+ from 0-5% most central collisions. As de- 
scribed in Eq. (3), the combinatorial background of the 
Q--K* correlation is modeled by the Q+-K+ correla- 
tion based on events with one Qt. Qt and K+ both 
containing 5 quarks and the narrow centrality bin make 
the Qt-K* correlation a good candidate for the combina- 
torial background in the Q--K* correlation. Similarly, 
we can also take the difference between Q+-K~ and Q7- 
K- to calculate the Q+-K~ correlation, which represents 
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FIG. 6. CBS correlations between Q7 (QT) and kaons in 0- 
5% most central Au + Au collisions at 7.7 GeV (|n| < 1). 
The left columns show the results of Cle a, and the right 
columns, those of Coe. 
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FIG. 7. CBS correlations between Q7(Q*) and kaons in 0- 
5% most central Au + Au collisions at 14.6 GeV. 


contributions only from SC. Therefore, any difference be- 
tween Q--K*+ and Q*+-K~ correlations will be sensitive 
to the BNT dynamics, which is absent in Scenario 2. 

Figures 6 and 7 show the CBS correlations, for (left) 
Q--K* and (right) Q+-K~ at 7.7 GeV and 14.6 GeV, 
respectively. The correlations are shown as a function of 
k*, rapidity difference (Ay), and transverse momentum 
difference (Apr) in three rows, respectively. At 7.7 GeV, 
the two AMPT versions exhibit clear differences. Com- 
pared with the SM version, the default version shows 
a stronger Q+-K~ correlation, presumably indicating a 
stronger and more localized SC. For the Q7-K* correla- 
tion, which is sensitive to both SC and BNT dynamics, 
while the total magnitudes (integral) are similar between 
the two AMPT versions, the correlation shapes are dif- 
ferent as a function of both k* and Apr. In the SM ver- 
sion, the coalescence formation mechanism and the BNT 
dynamics seem to yield a narrower correlation function 
between QT and K*. 
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FIG. 8. CBS correlations between Q7 (QT) and A(A) in 0-5% 
most central Au + Au collisions at 7.7 GeV (|n| < 1). The 
left columns show the results of 40 x CSES, and the right 


columns, those of OSPS. 


At 14.6 GeV (Fig. 7), the difference between the two 
hadronization schemes is relatively small. The shape dif- 
ference in the correlation as a function of k* and Apr 
is still visible, albeit much less prominent than at 7.7 
GeV. The CBS correlations at 7.7 GeV and 14.6 GeV sug- 
gest that the event-level Q7-K* correlation is stronger 
than the Q+-K~- one. Since + is only produced in Sce- 
nario 2, whereas Q~ can be produced in both scenarios, 
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FIG. 9. CBS correlations between Q7 (QT) and A(A) in 0-5% 
most central Au + Au collisions at 14.6 GeV (|n| < 1). The 
left columns show the results of 10 x CSES, and the right 
columns, those of OSPR, 


the stronger Q7-K* correlation indicates the presence of 
Scenario 1 for the Q~ production at both beam energies. 
Particularly at 7.7 GeV, Scenario 1 seems to contribute 
more to the Q7 production in the SM version, as the 
difference in correlation amplitudes between Q~-K* and 
Q+-K- becomes noticeably larger. 

Figures 8 and 9 show the CBS correlations between 2 
and A, and their anti-particle pairs as a function of three 
kinematic variables in 0-5% most central Au + Au col- 
lisions at 7.7 GeV and 14.6 GeV, respectively. Besides 
strangeness and baryon number conservation that govern 
the correlations between QT and A, the Q*-A correla- 
tions are also sensitive to BNT dynamics. At 7.7 GeV, 
the default version of AMPT displays stronger correla- 
tions than the SM version. 

For both AMPT versions, the Q~-A correlations are 
much weaker than the Q+-A ones, confirming that QF is 
only produced via Scenario 2. The default version also 
indicates a larger discrepancy between Q*+-A and Q7-A, 
which suggests a more prominent contribution from Sce- 
nario 2, complementary to the information on the differ- 
ence between the Q--K+ and Q+-K~ results. At 14.6 
GeV, the Q+-A correlations (the right column in Fig. 9) 
are also higher than Q~-A (the left column in Fig. 9) by 
about an order of magnitude. However, there seems to be 
no significant difference between the two AMPT versions 
at this energy. 
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0-5% most central Au + Au collisions at 7.7 GeV (|n| < 1). 
The left columns show the results of 5 x Oo ee. and the right 


columns, those of OSTE 


(a) } AMPT default 


“HMI seats : 
| A 


0.0 oS 1.0 IS 20 #25 3060 35 00 9s 10 I5 20 25 30 35 40 


k*(GeV/c) k*(GeVic) 
(c) (d) 


oo ih HE ead ae 


-2.0 -15 -10 -05 00 05 1.0 15 -2.0 -15 -10 -05 00 05 10 15 2.0 


Ay Ay 
0.015 (e) (f) 


m : 


4 tt 
p fame ynememseneag PR ni Hitters Epi ntnuntEnnnnnnnnnn 


—0.005 


0.0 0.5 1.0 hen, 2.0 2s ao a5 0.0 0.5 1.0 La 2.0 25 3.0 a5 4.0 
Apr(GeV/c) Ap7(GeV/c) 


FIG. 11. CBS correlations between Q7(Q+) and E*(E7) in 
0-5% most central Au + Au collisions at 14.6 GeV (|n| < 1). 
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columns, those of Co Fas 


Figures 10 and 11 show the CBS correlations between 
Q and © as a function of three kinematic variables in 0- 
5% most central Au + Au collisions at 7.7 GeV and 14.6 
GeV, respectively. At both beam energies, the stronger 
Q+-E- correlations relative to Q~-=+ also imply that 
Scenario 2 provides a larger contribution to the Qt pro- 
duction than to QT. At 7.7 GeV, the difference in the 
Q--5* correlations between the two AMPT versions sug- 
gests that the default version generates a larger contri- 
bution from Scenario 2 to the QT production. At 14.6 
GeV, the Q*+-=> correlations are still stronger than Q7- 
Z+ but no significant difference appears between the two 
versions. 


IV. SUMMARY 


The Q production in nuclear collisions at the RHIC 
BES involves dynamics of baryon number transport, 
strangeness conservation, and baryon number conser- 
vation. To investigate these effects, we have used the 
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